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1 Basic equations of magnetohydrodynamics 

A plasma is a mixture of ions, electrons and neutral particles. At the macroscopic level, beyond the kinetic 
description necessary in rarefied situations, dynamics of dense plasmas is governed by the interaction 
between the fluid components and the electromagnetic fields. When dealing with collisional plasmas [53] 
[57], it is often appropriate to use one-fluid approximations which are much simpler than a complete 
kinetic theory and relevant for a number of applications (astrophysics [20] [SO] |100j [104] , plasma physics 
m m, electrometallurgy [45] etc...). 

Moreover a further approximation is obtained when one studies the interaction between electromag¬ 
netic fields and a electric-conducting fluid with electrical effects negligeable compared to magnetic effects. 
This simplified scheme is called magnetohydrodynamics (MHD). It unifies classical fluid dynamics and 
magnetism of continuum media in a coupled way: electromagnetic fields induce currents in the moving 
fluid which produce forces which in turn modify the electromagnetic fields. If one considers a compressible 
heat-conducting fluid, the equations governing the system will be the compressible Navier-Stokes-Fourier 
system coupled to the set of Maxwell’s equations through suitable momenta and energy sources. 

1.1 Modelling of the magnetic field 

Recall that the electromagnetic field is governed by the Maxwell system [75] decomposed into 


the Ampere’s law 

dtT) + j = curhH, 

(1) 

the Coulomb’s law 

diVa;D = Pc, 

(2) 

the Faraday’s law 

-I- curRE = 0, 

(3) 

the Gauss’s law 

div,cB = 0, 

(4) 
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where D is the electric induction, B is the magnetic induction, E is the electric field, H is the magnetic 
field, Qc is the density of electric charges and j the electric current. We suppose that these quantities are 
linked by the constitutive relations 


B = and D = eE, (5) 

where the magnetic permeability and e the dielectric permittivity are positive quantities, generally 
depending on the electromagnetic field. Furthermore one observes that the electric density charge gc and 
the current density j are interrelated through the electric charge conservation 

dtQc + divj = 0, (6) 


however we assume in the following that electric neutrality holds: then g^ = 0. 

Finally we suppose that the magnetic induction vector B is related to the electric field E and the 
macroscopic fluid velocity u via Ohm’s law 


j = tT(E + u X B), 


(7) 


where the electrical conductivity a = a{g,'&,ll) of the fluid is a positive quantity. 

Now the magnetohydrodynamie approximation m [66] [98] consists in neglecting the term 9*0 in o 
and we obtain 

j = curl3,H. (8) 

This approximation corresponds to a strongly non-relativistic regime and has been justified in the 2- 
dimensional case by S. Kawashima and Y. Shizuta in |63j . and more recently in 3-dimensional space by 
S. Jiang and F. Li (see [57] for isentropic flows and [58] for the general case). 

Accordingly, equation ([3]) can be written in the form 


9tB -h curia; (B X u) -h curia. 


— curia; 

a 



= 0 , 


(9) 


where fi = ^(|H|) > 0. 
Setting 


equation ([3]) becomes 


7W(s) = f Tdr{Tg,{T)) dr, 
Jo 


at7W(|H|)+j-E = div, (HxE), 


( 10 ) 

( 11 ) 


Supposing now that we study the previous system in a bounded region 12 C we add boundary 
conditions 




( 12 ) 


and initial conditions 


B(0, x) = Bo(a:), E(0,x) = Eo(a;) for any x S 12. (13) 

In accordance with the basic principles of continuum mechanics, the magnetofluid dynamics is described 
(see Cabannes [12] 1 by the Navier-Stokes-Fourier system of equations with magnetic sources. This in¬ 
cludes mass conservation 

5te-kdiva;(£)u) = 0, (14) 

the momentum balance 


dt{gu) diva;(£)u 0 u) = diva;T -|- PcE -I-j x B -|- f, 


(15) 
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and the energy balance 

dt{^Q\vL\^ + ge'j + + ge +p)u- Su^ + div^jq = E • j + g, (16) 

where e is the specific internal energy, q denotes the heat flux, T is the Cauchy stress tensor, f is a given 
body force, 5 is a given energy source and the Lorentz force gcE+j x B is imposed by the electromagnetic 
field. In the magnetohydrodynamic approximation, it may be shown m that this approximation also 
implies that QcE, is negligible compared to j x B. Therefore, the momentum balance (flSl) reads 

9t(£iu) + diva;(£iu(g) u) = div 3 ,T + j x B + f. (17) 

Furthermore the Cauchy stress tensor T is given by Stokes ’ law 

T = S-pI, (18) 

where p is the pressure, and the symbol S stands for the viscous stress tensor, given by 

S = v (va;U + Vxu‘ - ^divx(u) 1^+77 diva;(u) I, (19) 

with the shear viscosity coefficient v, and the bulk viscosity coefficient rj. 

Finally we impose a no-slip boundary condition for the velocity 

Using CD and (fTCl) the energy of the system satisfies 

+ ge + 7W(|H|)^ -t div^, -t pe -h p) u-h E x B - Su q^ =0, (21) 

where e is the specific internal energy, and q denotes the heat flux. 

We prescribe a no-flux boundary condition for q 

q-n|an= 0 . ( 22 ) 

Using the previous boundary conditions we deduce that the (conserved) total energy of the system satisfies 

^ + pe - 1 - M(|H|)) da; = 0. (23) 

In the variational formulation it will be more convenient to replace (ED by the entropy balance 

dt{gs) + div 2 ;(psu) -t diva;(^^) = r, (24) 

with the specific entropy s, and the entropy production rate r, for which the second law of thermodynamics 
requires r > 0 . 

The specific entropy s, the specific internal energy e, and the pressure are interrelated through the 
second principle of thermodynamics 

■dDs = De + (25) 

where D stand for the differential with respect to the variable p, d, which implies some compatibility 

conditions between e and p ( Maxwell’s relation). 

If the motion is smooth, the sum of the kinetic and magnetic energy satisfies 

dt -I-A4(|H|)^ -I- div,;, -|-p^u-|-ExB - Su^ = p divx(u) - § : Va,u - ^|curla;B|^. 

(26) 
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Indeed using (1^ one checks that 


q- 

■d 


(27) 


1 


1, 


r = - ( S : Va:U+-|curl:j;H| - 


If the motion is not (known to be) smooth, validity of (fTOl) is no longer guaranteed and the dissipation 
rate of mechanical energy may exceed the value § : Va;U. 

Accordingly, we replace equation (l24l) by inequality 

dtigs) +div,c(psu) + div^; + i|curla;Hp - , (28) 

to be satisfied together with the total energy balance (l23ll . 


1.2 Hypotheses on constitutive relations 

In order to close the system one needs first to postulate a constitutive equation relating the pressure 
to the other variables. Several possibilities have been considered. We choose here, in the spirit of [50] 
(however see [36] [39] [51] for other possible choices) to focus on high-temperature phenomena appearing 
in various astrophysical contexts (stars, planetar magnetospheres, neighborhood of black holes etc...) [00] 
m CHD] cni, or plasma physics applications (tokamaks, inertial confinement fusion devices etc...) [64] 

EZ]. 

In that case, the so called equilibrium diffusion limit of radiation hydrodynamics involves a modified 
pressure law where the pressure in the fluid is augmented by a radiation component pr{'&) related 

to the absolute temperature through the Stefan-Boltzmann law 


Pr{'&) = with a constant a > 0 


(29) 


(see, for instance. Chapter 15 in [33]). Accordingly, the specific internal energy of the fluid must be 
supplemented by a term 

enW = eRiQ,^) = (30) 

9 

where g is the fluid density and similarly, the heat conductivity of the fluid is enhanced by a radiation 
component (see [IT] [75]) 


qi?(^^) = —with a constant kr > 0 (31) 

Finally the total pressure p takes the form 

p{g,'d)=pF{g,'d)+PR{'d), (32) 

where the radiation component is given by (HU). Similarly, the specific internal energy reads 

e(£»,'d) = CF (£»,■(?) +e/i (£»,■(?), (33) 


with cr determined by dSU. 

Furthermore we suppose that the gas is perfect and monoatomic 

2 

PF{g,'&) = ■^geF{g,'&)- 
The following hypothesis expresses the convexity of free energy 


dpF{g,'&) 

dg 


deF{g,'&) 

dd 


> 0 for all p, -d > 0, 


(34) 


(35) 
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and finally, we suppose 


liminf 1 ?) = 0, liminf > 0 for any fixed > 0, 

s^o+ ^ ^ g^Q+ dg 


and 


liminf ep(g, d) = 0, liminf > q f-Qj. gxed g > 0. 

t?->-o+ iS-^os- og 


As one can check by solving (l25l) with (l34l) . there exists a function Fp G C^(0,oo) such that 

PF(g,i 9 ) = i 9 ^Fp(^). 

Consequently 


so according to (I55|) 
where the latter inequality yields 

In particular 


Pp{z) > 0, -Pp{z) — Pp{z)z > 0 for any z > 0, 
o 


(W)' 

V ZS / 

Pp{z) 


< 0 for all z > 0. 
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Poo > 0 for z —>■ oo. 


According to (l36t . (l37)) . we see that Pp G C'^[0,oo). Moreover 


0 < liminf - (-Pf(z) — PP(z)z) < limsup - f-PF(-2:) —-Pf(-z)2:') < oo, 

z—¥Q-\- 2 V u / 2 —>^ 0 + 2 ^ ^ O / 

1/5 \ 

limsup-(-P f( 2:) - Pf < oo, 

5::—VriO 2^0 / 


and 


^->■o+ 

Now it follows from (j^H) that 


P' (z) 5 

lim Pp{z) = 0, lim Pp{z) > 0, lim —q— = -poo > 0. 

■ " ' . n , z— >00 2 ; 3 3 


z-s-OG 


Ppiz) > CZ3 for all z > 0, and a certain c > 0. 
Therefore there exists Pc > 0 such that the mapping 

g^ PF{g,'&) -Pcg^, 

is a non-decreasing function of g for any fixed d > 0. 

Accordingly, using equation (1251) we have 


s{g, '&) = SF{g, 1 ?) -f SR{g, -d), 


where 


and 


SF{g,d) = Sf(-^), with S'^(z) 

\li2 / 


3 iPpiz) - P'f{z)z 


QSR{g,'&) = 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


(45) 

(46) 

(47) 

(48) 
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Note that after ((5^ and (li^ . Sp is a decreasing function such that 

lim = -^-Pf(O) < 0, 

z->'0+ 5 

whence, normalizing by 5'ir(l) = 0, we get 

— Cl log(z) < Spiz) < —C 2 log(z), Cl > 0, for all 0 < z < 1, (49) 

and 

0 > Sf{z) > —C 3 log(z) for all z > 1. (50) 

Using the second law of thermodynamics, we see that viscosity coefficients iz and rj are non-negative 
quantities. 

In addition, we shall suppose that 

i. = |H|) > 0, 77 = r 7 (ft, |H|) > 0, (51) 

where the bulk viscosity v satisfies some technical, but physically relevant, coercivity conditions to be 
specified below. 

Similarly, the heat flux q obeys Fourier’s law 

q = qi? + qF, (52) 

where the radiation heat flux is given by (1311) and 

qF = -7t;^(e,ft,|H|)V,d, F>0. (53) 

For the sake of simplicity we assume the transport coefficients kf, and a to admit a common 
temperature scaling, namely 

ci(l + ft“)<FF( 0 ,i?,|H|), Kft,|H|),a-i( 77 ,ft,|H|) <C 2 (l + ft“), Cl, C 2 > 0, (54) 

with a > 1 to be specified below. Similarly, we suppose that 

0<77(ft,|H|)<C3(l + 7?“), (55) 

and that the magnetic permeability has the symbol property 

Cj,s(l -I- s)~'^ < dg{sfj.{s)) < Cfcs(l -I- s)~^, (56) 

for any s > 0 and for fc = 0,1, with Ck > 0. 

We finally suppose that external sources are absent: f = 0 and g = 0. 

Let us conclude this introduction by briefly fixing the limits of our study. As our framework is 
explicitely viscous and compressible, all the important works about the incompressible model are beyond 
our scope. However the reader may consult, among standard references: Ladyzhenskaya and Solonnikov 
[69], Duvaut and Lions m m, Sermange and Temam [88] and Cannone et al. m for reviews of 
the incompressible setting and Germain and Masmoudi [46] (with references therein) for the inviscid 
case. Finally we refer the reader to Eringen and Maugin [34] for more complex situations coupling 
electrodynamics to continua. 

The article is structured as follows: in Section |2| we consider the main existence result of weak 
solutions, in Section|3|we present related models (barotropic or polytropic) and in SectionUjwe introduce 
asymptotic models. In Section [5] we present the main results concerning existence and uniqueness of 
strong solutions including critical aspects and finally results for one dimensional models are given in the 
last Section m 
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Weak solutions 


2.1 Variational formulation 


As in the pure hydrodynamical model [52] the approach is based on the concept of variational solutions. 

1- The mass conservation is expressed through the weak form 





+ gu ■ 


dx dt + 


go(p(0,-) dx = 0 , 


(57) 


to be satished for any test function p G 27([0,T) x R^) where go is the initial density. 

Assuming that both g and the flux giu are locally integrable on [0,T) x 12 equation (l57ll makes sense. 
Moreover if g belongs to L^(0,T; L^(fi)) for a 7 > 1, then [82] g € C([0,T]; i.e. 


g(t) —)> go weakly in L^(fl) as t —>■ 0 , 


and the total mass 


Mo= g{t) dx = / £>0 da;, (58) 

Jq Jq 

is a constant of motion. 

For 7 > 2 one can use the DiPerna-Lions theory [55] to show that the pair (g, u) satisfies the renor¬ 
malized equation 



(b{g)dtp + b{g)u ■ V+ {b{g) 


b'{g)g)diya: 



dx dt -I- / b{go)p{0^ •) dx = 0 , 
an 


(59) 


for any p G X>([0,T) x R.^), and for any continuously differentiable function b whose derivative van¬ 
ishes for large arguments. It can also be shown that any renormalized solution g belongs to the space 
C'([ 0 ,T];Li(^^)) (see [28]). 

2- The variational formulation of the momentum equation reads now 


J j ^(pu) • dtp + (eu ® u) : y^^p+p diVa, (</;)) dx dt 

[ / - (j X B) • (£)') dx dt - [ (£iu)o • p{0, ■) dx, 

Jo Jn ^ ' an 


(60) 


which must be satisfied for any vector field p G 27([0,T) x 12; R^). Once more we suppose that the 
quantities pu, £iu ® u, p, S, gVxdf, and J x B are at least locally integrable on the set [0, T) x 12. 

The initial value of the momentum (£<u)o must satisfy 


(£)u)o = 0 a.e. on the set {£>0 = 0 }. 
Finally using the no-slip boundary condition (1501) we get 

uG L2(0,r;W'o’^(12;R3)). 


(61) 

(62) 


Observe that our choice of the function space has been motivated by the a priori bounds resulting from 
the entropy balance specified below. 

3- According to (155)) the weak form of the entropy production reads 



(^gs dtp 


+ gsu ■ WxP + ^ ■ '^xP 


dx dt < 


(63) 


f 1 /q- 

io in ^ V 


— § : VxU-^|curl 2 ;Bp ')p dx dt 

J 


(£is)o;£>( 0 , •) dx, 
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for any test function ip S X>([0,T) x R^), yi > 0. 

We observe that the presence of d in the denominator forces this quantity to be positive on a set of 
full measure for (l63l) to make sense. 

4- Finally Maxwell’s equation Q is replaced by 

J j X u + — curia; —^ • curla;93^ da: dt + J Bq • i^(0, •) da;, (64) 

to be satisfied for any vector field ip £ 2?([0,T) x R.^). 

Hereafter using boundary conditions m and (15^ one must take 

Bo £ diva; (Bo) = 0 in Bq • n|an = 0. (65) 

By virtue of Theorem 1.4 in [^, Bq belongs to the closure of solenoidal fields on 'D{il) with respect to 
the norm. 

Applying (l23l) and (l28l) . we infer that A1(|H|) £ L°°(0, T; L^(H; R^)) andcurla;H £ L^(0, T; L^(H; R^)). 
But using (fTOl) and (IMl) we deduce that B £ L°°(0, T; L^(r2; R^)), and from (jU and (IT^ we see that 

diva;(B)(t) = 0 in ^'(B), B(t) • n|ao = 0 for a.e. t £ (0,T), 

so using Theorem 6.1 in [32] we conclude that 

Hg diva;(B)(t) = 0, B • n|an = 0 for a.e. t £ (0, T). (66) 

5- According to (1311) we shall assume the total energy to be a constant of motion 

E{t)= J (^ip(t)|u(t)p + (£ie)(t)+7W(|H(t)|)j dx = Ao for a.e. t £ (0,r), (67) 

where 

Ao = y (^^|(£<u)oP + (ee)o + A4(|Ho|)) dx. (68) 

Of course the initial values of the entropy {gs)o and the internal energy {ge)o should be chosen consistently 
with the constitutive relations. A natural possibility consists in fixing the initial temperature i9o and 
setting 

(gs)o = gos(go,^o), (ge)o = goe(go,i^o)- (69) 

2.2 Global existence of weak solutions for large data 

The following result holds 


Theorem 2.1 Let B C R^ be a bounded domain with boundary of class 5 > 0. 

Suppose that the thermodynamic functions 

p = p{g,-d), e = e{g,-d), s = si^g^d) 

are interrelated through (I2S1), where p, e can be decomposed as in dSlI), (|33|), with the components pp, 
ep satisfying dMj). Moreover, let pp{g,'t)), ep{g,'d) be continuously differentiable functions of positive 
arguments g, tt satisfying (|351)-(|37|1. 

Furthermore, we suppose that the transport coefficients 


V = v(ff, |H|), T] = piff, |H|), Kp = Kp{g,-d, |H|), cr = a{g,d, |H|) and p = Ai(|H|), 


are continuously differentiable functions of their arguments obeying with 


1 < a < 


65 


Finally, let the initial data qq, (^?u)o, i9o, Bq be given so that 

go e L^n), {gu)o £ L\n-,R^), do G Bq G L\n-,R^), 


(70) 


(71) 


^>0 >0, do > 0, (72) 

{gs)o = gos{go,'do), —|( 6 iu)o|^, {ge)o = goe{go,do) e (73) 

Qo 


and 

div 2 ;(Bo) = 0 in V'{n), Bq • n|an = 0. (74) 

Then problem (I571) - (I55]) possesses at least one variational solution g, u, d, B in the sense of Section \2.1\ 
on an arbitrary time interval {0,T). 


Observe that hypothesis is a technical restriction. However functional dependence on d was rigorously 
justified in [26] after an asymptotic analysis of some kinetic-fluid models. The same temperature scaling 
on all transport coefficients was also suggested in [47]. 

2.2.1 Difficulties and methods 

The main difficulty of the problem being the lack of a priori estimates, methods of weak convergence based 
on the theory of the parametrized (Young) measures represent the main tool [83] to proceed. Accordingly, 
we denote by b{z) a weak limit of any sequence of composed functions { 6 (z„)}„>i. More precisely, 

/ b{y,^)T{y)<iy= ‘p{y)( b{y,z) dAy{z)) dy, 

where Ay{z) is a parametrized measure associated to a sequence {z„}„>i of vector-valued functions, 
z„ : R^ —>■ (see Chapter 1 in [83]). 

Assume there is a sequence of approximate solutions resulting from a suitable regularization process. 
The starting point is to establish the relation 

(^p{g,d) - (^^:/-k 7 ?)div,j,(u)^ b{g) = (^p{g,d) b{g) - |B|) -f 77 ( 1 ?, |B|))diV:j,(u)^ b{g), (75) 

for any bounded function b. 

The quantity p— ( 4 / 3 jz-|- 77 )diva;(u) is called the effective viscous pressure and relation (IT^ was proved 
by P.-L. Lions m for the barotropic Navier-Stokes system, with p = p{g), when viscosity coefficients z/ 
and r] are constant. The same result for general temperature dependent viscosity coefficients was obtained 
in [40] with the help of certain commutator estimates in the spirit of [T^ . 

In the present setting, this approach has to be modified in order to accommodate the dependence of 
iz and r] on the magnetic field B. 

The propagation of density oscillations can be described by the renormalized continuity equation [28] 
dtb{g) + diva,( 6 (p)u) -k (b'{g)g - b{g)'j diva;(u) = 0, (76) 

and its “weak” counterpart 

dtb{g) + dwa:{b{g)u) + (b'{g)g - b{g)^ div,j,(u) = 0. (77) 
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To prove ( 1751 ) one must show first boundedness of the oscillations defect measure 

osc^+i[p„ -)> q]{Q) = sup (limsup [ \Tk{Qn) - T’fc(e)p+^ da; dt) < oo, (78) 

fc>l n^oo JQ ' 

where Tfc(p) = minjp, k}, for any bounded Q C K'* and a 7 > 1. 

Due to the poor estimates resulting from (l54l) . (l55l) and (l63l) relation (l78l) has to be replaced by a 
weaker weighted estimate 

sup (limsup / (1 + 'dn)~^\Tk{gn) - da; dt) < 00, (79) 

fc>l n—^oo JQ ’ 

for suitable ^ > 0, 7 > 1. 

Finally one recovers strong convergence of the sequence {i?ri}n>i of approximate temperatures, know¬ 
ing that spatial gradients of dn are uniformly square integrable and that 

gs{g,'d)'d = gs{g,'&)'d, (80) 

where ( 1551 ) can be deduced from the entropy inequality (1631) by using Aubin-Lions lemma. 

The structure of the proof of Theorem 12.11 goes as follows: introducing a three level approximation 
scheme adapted from [35] allows to reduce the proof to a weak stability problem. One first proves 
uniform bounds on the sequence of approximate solutions. Relying on these uniform bounds, one can 
handle the convective terms in the field equations by means of a compactness argument. The strong 
convergence of the sequence of approximate temperatures is proved by using entropy inequality together 
with the theory of parametrized (Young) measures discussed below. Then one finally shows that the 
sequence of approximate densities converges strongly in the Lebesgue space T^((0, T) x D). This requires 
weighted estimates of the oscillations defect measure in order to show that the limit densities satisfy the 
renormalized continuity equation. 

2.2.2 The approximation scheme 

The approximation scheme used in the following is inspired from Chapter 3 in [42] . 

1- The continuity equation dm is regularized by a viscous approximation 

dtg + div2;(£(u) = eAg, e > 0, (81) 

on (0,T) X n, supplemented by the homogeneous Neumann boundary conditions 


The initial density is given by 
where 


n|an = 0. 


ll(0; *) P0,(5, 


(82) 

(83) 


go,s G (7^(11), Vxgo,s ■ n|an = 0, inf go,six) > 0, (84) 

with a positive parameter 5 > 0. 

The functions go,s are chosen in such a way that 

go,5 go in At(D), |{po.5 < eo}| 0 for (5 0. (85) 

Here, of course, the choice of the “critical” exponent 7 = 5/3 is directly related to estimate (1331) established 
above. 
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2- The regularized momentum equation is 

dt{gu) + div 2 ;(eu(g) u) + VxP + +eVa;uVa;£» = div 3 ;(S) + j x B, ( 86 ) 

in (0,T) X n where the approximate velocity field satisfies the no slip boundary condition 

u|aa = 0. (87) 

We also prescribe initial conditions 


where 


(£»u)o,5 = 


( 0 U)(O, •) = (£»u)o, 5 , 

(£)u)o provided go,s > go, 
0 otherwise. 


( 88 ) 


(89) 


The artificial pressure term in () 86 l) provides additional estimates on the approximate densities in 
order to facilitate the limit passage £ —>■ 0. To this end, one has to take T large enough, say, T > 8 , and 
to re-parametrize the initial distribution of the approximate densities so that 


S f £Iq ^ dx —>■ 0 for ^ 0. 

Jn 

3- The entropy equation IHM is replaced by the regularized internal energy balance 
dt{ge + 6g'd) + diva, (^{ge + Sg'd)u^ — div^, (^{kf + 


1 


2 „r -2 


= S : Va;U - p diva;(u) + - Icurl^-HI + er\V,^g\^g- 

a 

to be satisfied in (0,T) x fl, together with no-flux boundary conditions 

Vj;!? • n|a = 0. 

The initial conditions read 

g{e + 5-d)(0, •) = go.sie{go,s, '&o,s) + S-do^s), 
where the approximate temperature satisfies 

'do,s G ■ n|an = 0, inf i?o.5(a;) > 0, 




and 


1 ^ 0 ,s —t 'do in LP{n) for any p > 1 , 
d go,s log(i9o,s) dx -)> 0 , 

Jq go,ss(go,s, 1 ^ 0 , 6 ) dx ^ gos(go,^o) dx, 


as <5 —)> 0 , with 


/ go,se(go,s, 1 ^ 0 , 6 ) dx < c uniformly for d > 0 . 

Finally, the magnetic induction vector B obeys the original Maxwell’s equations 

dtB + curla;(B X u) + curl^, [ —curl^, ( — ) ) = 0, diva;B = 0 


(90) 


(91) 


(92) 


(93) 


(94) 


(95) 


(96) 


(97) 
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in (0,T) X Q, supplemented with the initial condition 


B(0,-) = Bo.5, (98) 

where, by virtue of Theorem 1.4 in |93] . one can take 

Bo,( 5 G diva;Bo,5 = 0, (99) 

Bo ,5 ^ Bo in ^ ^ 

For given positive parameters e, S , and F > 8, the proof of Theorem 12.II consists in the following steps: 
Step 1 Solving problem (I5T]) - (I100I) for fixed e > 0, <5 > 0. 

Step 2 Passing to the limit for e — >■ 0. 

Step 3 Letting i5 —>■ 0. 


One sees that Step 1 can be achieved by using a simple fixed point argument: equation (IM)) is solved 
in terms of the velocity u with the help of the Faedo-Galerkin method, where g, 'd, and B are computed 
successively from (EH), EID, and (EZl) as functions of u. It is easy to check that the corresponding 
approximate solutions satisfy the energy balance 


d 

dt 


L 


pe + 7W(|H|). 


F- 1 


+ Sgdj 


dx = 0, 


in T>'(0,T), and the limit property 


‘^Qo,s 


lim f f-£»|u 
i->0+ Jq\2 

I(£’u)o.5 


■pe + 7W(|H|) 


r -1 

^ + Qo,se-igo,s, 'do,a) + -^(1110,51) 4 


g^ + Sg'd^ 

6 p 

P _ ^ 00,s 


dx = 


( 101 ) 


( 102 ) 


^9o,s'&o,s^ 


dx. 


The reason for splitting Steps 2 and 3 into the e and (5—parts is a rehned density estimates based on 
multipliers Va;(—A)“^[p^], where one must take ,5 = 1 when the artificial viscosity term eAg is present, 
while uniform estimates require /3 to be a small positive number. For this reason, we focus only on Step 
3: the goal is to establish the weak sequential stability property for the solution set of the approximate 
problem: 

1. The density > 0 and the velocity satisfy 


/ / (gsdtip + gsus-'Vx‘p)dxdt+ / po,5<p(0, •) dx = 0 

Jo Jn ^ 

for any test function Lp € 22([0,T) x R^). In addition, 

us(t) G VFQ’^(fl;R^) for a.e. t G (0,r). 


(103) 


(104) 


2. The momentum equation 

i-T 


P5U5 • dtip 4 (p5U5 (g) Ui) : 4p5 diva;((p) 4 Sgg div^ 


dx dt = 


(105) 


Jo Jfi ^ 



Ss ■■ Vx(p 


(ji X Bi) -(f^ dx dt- J (pu)o,5 • ip{0, ■) dx 


holds for any ip G T’([0,T) x n;R^), where ps = p(p 5 ,d 5 ), and S^, are determined in terms of u^, 
and B^ through the constitutive relations. 
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3. The entropy production inequality reads now 



(^{gsss + Sgs log(i9s))dtip + (gsssus + 6 gs \og{ds)us) ■ dx dt 


(106) 


+ ^ ^ - r5<-V.da • dxdt 

< [ - §5 : Va;U5 - — |curla;H5|^ - dx dt 

Jo Jn vs\ ds as ° J 


■ J (^Qo,6s{go,s, '&o,s) + 6go^s log(do,5)) (/?(0, •) dx, 


for any (p G T>([0,T) x R^), p > 0. Here, da is assumed to be positive a.e. on the set (0,T) x H, 
S (5 = s{gs,'ds), as = a{gs,'ds), and is a function of gs, tJs given by (|5^ . 

4. The magnetic induction vector satisfies 


/ / (b^ • i9t(p - (B^ X Ui + — curl2,H5) • curUv?) dx dt + / Bo,5 • (/^(O, •) dx = 0, (107) 

Jo Jn ^ ' Jna 


for any vector field p G 27([0, T) x K^; 

5. Finally, the (total) energy equality 


J + gse + JAiHs) + + Sgs'Js'jit) dx 


(108) 


/ 1 


\‘2go,( 


■|(eu)o,5p + go,se{go,s,'&o,5) + A4(Ho,5) dx 


'^6 


r-1 


£'o 5 + ^6'0.5^?0.5 dx 


holds for a.e. t G (0,T). 

The weak sequential stability problem to be addressed below consists now in showing that one can 
pass to the limit 

gs -t p, Us ^ u, -ds J, B^ B, 


as (5 —>■ 0, in a suitable topology, where the limit quantity {p, Ujd, B} is a variational solution of problem 
(|57|l-l6l 


2.2.3 Uniform bounds 

One must identify uniform bounds on the sequences {P( 5 }, 5 > 0 ’ {'^<5}(5>0’ {'d< 5 } 5 >oi ^^d {B^j^^Q through 
total energy balance pOSp . dissipation inequality p06p and the other relations resulting from p03p - p08p . 

1. Total mass conservation. As gs, gsus are locally integrable in [0,T) x O, it follows from (11031) that 
total mass is a constant of motion, specifically, 


/ gs{t) dx = / go^s dx 


Jn Jn 

In particular, as > 0 and (|^ holds, we get 

{£'<5}(5 >o bounded in 

2. Energy estimates. Using (IMl) and (l44ll one gets 


: Mq for a.e. t G (0, T). 

L°^{0,T;L\n)). 


gse{gs,iJs) dx = - / PFigs,'ds) dx > 


3po 


Qs da:- 


(109) 


(110) 
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Total energy balance (11081) together with bounds on the initial data (1551) . (IMl) . (TOl) . (IMl) . and (llOOl) yield 
the energy estimates 

{e5}5>o bounded in L“(0,T;L3(r2)), (HI) 

bounded in L°°(0,T;L'^(0)), (112) 

(113) 

(114) 


bounded in L“(0,r;L2(f7;] 

, {£'5e(£>5,i?5)}5>0 bounded in L°“(0,T;L^(n)), 


and 


6 dec < c uniformly with respect to d > 0. 

Jn 

In particular, using Holder’s inequality, (Hill) and H14I) imply 

{(pu) 5 }^^p bounded in L°“(0, T; (H; K^)). 


(115) 


(116) 


3. Dissipation estimates. Taking a spatially homogeneous test function ip such that p{0, •) = 1 in the 
entropy production inequality H06I) and using (IMl) (IMl) . we obtain 




V a: 1^(5 


q<5 ■ Vaj'di , 1 


■ds 


+ —Icurl.Hap + dx dt (117) 

as J 


< Cl + / £' 5 (t)s(p 5 (t),i 95 (t)) + dp 5 (T) log(d 5 (T)) dec < C2 + / P5 (t)s(p5(t),i?5(t)) dx, 

JO J Vt 

for a.e. r S ( 0 ,T), where the last inequality follows from (jllip . pi 2 p . 

Using (| 4 ^ and ([ 50 p . the integral in the right hand side may be estimated 

[ gss{gs,'ds) dx = [ + gsSp ( -^ ) dx < ci - C2 / ^ gs (log(£»5) - | log(r?5) ) dx 

JQ. Jnd K-dl J V 2 y 

<03 + 04/ 3 p4log('d4) dx < C5 + C6 / gs'0sdx<C7, ( 118 ) 

Jtps<^?} -Jn 


where we have taken into account Hill) and H12I) . 

Thus the integral in the left hand side of H17p is bounded independently of 6 and we conclude after 
(I55P - (I5SP that 

|(l + d 5 )“ 2 “ (Va;U 4 )|^ is bounded in L^(0, T; L^(r2; K^’^^)), 


|(l + i?5) 2 divj; (u5)|^ is bounded in L^(0,T; L^(r2)), 

{V3; log(i?5)}5>o , are bounded in L^(0,r;L^(H;R^)), 

L J ( 5>0 

|(l + i?5)“2-curUH^I is bounded in L^(0,T;L^(U;]R^)), 


s>o 


and 


is bounded in T^(0, T ; T^(H; M^)), 
where a satisfies dTO]), and where we denote by 


(119) 

( 120 ) 
( 121 ) 
( 122 ) 

(123) 


") - -trace(Q)I, 


the traceless component of the symmetric part of a tensor Q. 
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Since the velocity field satisfies (11041) we find 


bounded in L^(0,T;WQ^’^(n;]R^)), (124) 

and (I121D and (I112p lead to 

|i9|| bounded in (^25) 

I J (5>0 

4- Positivity of the absolute temperature. According to physics, temperature must be positive a.e. on 
(0,T) X n. To justify that we use the uniform L^— estimates of Va;log(-d5) given by (I124|) and the 
following version of Poincare’s inequality: 


Lemma 2.1 Let LI C be a bounded Lipschitz domain, and uj > 1 be a given constant. Furthermore, 
assume O G LI is a measurable set such that |0| > o > 0. Then 

Ikllwu2(n) < c(o, fl,cc)(||V xv\\L^{n.R3) + da;)‘^y 


In order to apply Lemma l2.11 we show first the uniform bound 

0<T< / dx for a.e. t G (0,T). 

Jn 

As a consequence of (110611 we get 

J (^gss{gs,^s) + Sgslogiidsi^it) ^ J (^Qo,6s{go,s,'do,s) + 5go,s^og{'do,6)'^ dx, 

where, by virtue of (IM|) and (IMll 

6 / Po.5 log(i5o,5) dx0 and / go,s s{go,s,'do,s) dx ^ / £>os(e'o,'do) dx. 
On the other hand (lllip and ()112l) yield 

S / log('d5)(t) dx < 5 / gs{t)'ds{t) dx ^ 0 for a.e t € (0, T), 

Jq Jq 

while, using (|46l) 


J gssigs,'0s){t) dx = J ^ 


Supposing the opposite of (|126L one could extract a sequence l^at 

dsits) —>■ 0 weakly in L^(f2), and strongly in LP{LI) for any 1 < p < 4, 


where 


( + 9oSf ( I I dx < liminf [ gs{ts)SF(^)(ts) dx, 

gs{ts) ^ g{t) weakly in L^Ll), / gs{ts) dx = / dx. 

Jq, Jq 


(126) 


(127) 

(128) 
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Now, for any fixed K > 1, one can write 


gs{ts)SF(-^;) (ts) dx 

\ .Q 2 / 


(129) 


< 


' {Q5{ts)<K'Sj {ts)} 
< C I 


gs(ts)SF ( ) i^s) dx + f 3 g5{ts)SF ( 

\i9| / ''{e5(ts)>Ki}](ts)} \i?| 

£ 15 (^ 5 ) 11 + ^(4) j dx + SFiK) j go,s dx 


{ts) dx 


'{g5(ts)<K^^ (tj)} 


< 2c / dg {ts) dx + Sf{K) / go ,<5 dx. 

Jn JQ 

Thus combining (I127I) - (I129I) one concludes 

/ + QoSf(-^^ dx < Sf{K) [ go dx for any K > 1, 

dn 3 ^1^0 

which is in contradiction with the fact that Sf is strictly decreasing with limi<-_,.oo Sf{K) < 0, and 1)0 
positive (non-zero) on f2. This proves the uniform bound (I126I) . 

Finally, observing that 

T-e|n|< [ i? 5 (t) dx < |{i? 5 (t) > £}|^||i?5(t)||i,4(n) for any e > 0, 

J{-ds{t)>e} 

we deduce, using (I112L that there exist e > 0 and o such that 


|{i? 5 (t) > ell > o > 0 for a.e. t e (0 ,T), 

uniformly for (5 > 0. Then using (11211) one can apply Lemma |2.II to log(i? 5 ) and we obtain the required 
estimate 

{log('d 5 )}^^Q bounded in L^(0,T;VF^’^(n)). (130) 

5. Estimates of the magnetie field. Using the integral identity (I107|) we get 

div,(B 5 )(<) = 0 in V'{n), Ba • n|aa = 0, (131) 

which, together with estimates (11131) . (11221) and Theorem 6.1 in [32], yields 

bounded in L^(0,T;1U^’^(U;R^)). (132) 

6 . Pressure estimates. We first observe that estimates (llll|) . (11241) imply boundedness of the sequences 

{giU^I^^P and {g 5 U 5 0 in ^^*((0, T) x Vt) for a certain p > 1. 

Now one gets 


(VxUj) + \/'dsr]{'ds)^divx{us), (133) 

where, by virtue of (I112|) . (11191) . (11251) . and hypothesis (iTOl) . the expression on the right-hand side is 
bounded in L^((0,T) x 11) for a certain p > 1. 

Finally the Lorentz force x B^ = curla;H 5 x /^(H^)!!^ is controlled by (11131) . (I132p and (1531) which, 
combined with a simple interpolation argument, yields 

{ji X Bij^^g bounded in L^{{0,T) x 12) for a certain p > 1. (134) 
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At this stage, using the main result of [U], we are allowed to use the quantities 

= 'tlj{t)B[g‘g], iIj G 'D{0,T) for a sufficiently small parameter w > 0, 
as test functions in (jl05|l . where B is the Bogovskii operator [9] solving the boundary value problem 


divj; 





B\dn = 0 . 


(135) 


The resulting estimate reads 

fT 


'0 


j {p{QS)'^s)dt + dx dt < c, with c independent of S, 


then 

and 


{p{gs, ^^( 5 )} 5 >o is bounded in L^{{0, T) x fl) for a certain p > 1, 
is bounded in L^((0,T) X n). 

^ J ( 5>0 


(136) 

(137) 

(138) 


2.2.4 Sequential stability of the equations 

1. Continuity equation. With the estimates established in the previous section, we can pass to the limit 
(5 ^ 0 in (fnm . 

Indeed, extracting subsequences if necessary, we deduce from (I111|) . (I116|) . (11241) . and the fact that gs 
and U 5 satisfy the identity (110311 

Pi^pinC([0,T];LL,(ll)), (139) 

U 5 — ^ u weakly in L^(0, T; WQ^’^(n; M^)), (140) 

PiUi —^ pu weakly-(*) in L°°{0, T ; (11; R^)), (141) 

where the limit quantities satisfy dSTl). 

2.Momentum equation. Using the estimates obtained in Section [2.2.31 7 together with (I105|) we have 

gsns ^ gu in C{[0,T]-, ci (142) 

P 5 U 5 (g) u^ —>• gu (g) u weakly in L^(0, T; (11; R^^^)), (143) 

where we have used the embedding Wg^’^(H) ^ L®(H). 

Furthermore, in accordance with (I136|) . (11371) 

p{gs, -ds) pig, ’d) weakly in T^((0, T) x H), (144) 

and 

Sg^ ^ 0 in ^^’((0, T) x H), (145) 

for a certain p > 1. Here, in agreement with Section r2. 2. 31 

[ [ pig,'d)(p dx dt = [ / <pf / p(p,-d) dAt,a;(p,-d)') da; di, e 21((0,T) x H), 

Jo Ja Jo Jn ^ 

where At^xig, '&) is a parametrized (Young) measure associated to the (vector valued) sequence {gs, 
Similarly, one can use (110711 together with estimates (11131) . (11321) to deduce 

Hi ^ H weakly in 7.2(0, T; 1U^’^(11;R^)) and strongly in L^ii0,T) x 11;R3), (146) 
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and 

^(Hi) —>• /^(H) strongly in L^((0,T) x r2;R^), 


then 

for a certain p > 1. 


X Bj —>■ j X B weakly in L^((0,T) x r2;K^), 

Thus the limit quantities satisfy an averaged momentum equation 

J j (^(pu) • dtp + {gu (g) u) : + p{q, ■&) divx(</ 5 )^ dx dt 


(147) 


= J y (s : - (j X B) • dx dt- J {gu)o ■ p{d, ■) dx 

for any vector field p S 27([0,T) x n;R^). 

The symbol S denotes a weak limit in T^(0, T ; R^y^)), p > 1 of the approximate viscosity tensors 

Sa specified in (I133p . Clearly, relation (11471) will coincide with the (variational) momentum equation (1501) 
as soon as we show strong (pointwise) convergence of the sequences {£' 5 }, 5 >o {^< 5 } 5 >o- 


2.2.5 Entropy inequality and strong convergence of the temperature 

1. Entropy inequality. In order to obtain information on the time oscillations of the sequence 
we use estimates on dt{gss{QsT'ds)) provided by the approximate entropy balance (11061) . 

It first follows from (H5)) . (1501) that 

|£< 5 S( 05 , 1 ? 5 )| < + gs\ log(£> 5 )| + gs\ log(i? 5 )|), 

therefore, by virtue of the uniform estimates (IlllL (I112L (11161) . and (11301) . we can assume 


gss{gs,'&s) -t gsig,^) weakly in LP(( 0 ,T) x n), 
gss{gs, ^ 5 ^)u 5 — 7 > gs{g, i?)u weakly in T^((0, T) x Q-, R^), 


for a certain p > 1 . 

Similarly, one estimates the entropy flux 


^ < c(|V, \og{^s)\ + , 


where 


2 










«)2_ 


Choosing the parameter 0 < s < 1 small enough so that 


( 1-4 = 1 , 


we have after Holder’s inequality 

II ■■ 


LP(n;R3) — 




L2(n;B3) 


L4(n) 




Le(n) 


Then using estimates (11121) . (11231) together with the imbedding IT^’^(H) ^ T®(11), we get 

Vx'di —^ 0 in L^((0,T) x r2;R^) for a certain p > 1. 


(148) 

(149) 


(150) 


(151) 
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Moreover using similar arguments one can also show that 


I is bounded in L^((0, T) x R^), for a certain p > 1. 

J 5>o 


On the other hand according to (11211) 


b{'&s) Kd) weakly in L‘^((0,T) x ft)), and weakly in L^(0,T; kb^’^(r2)), 

for any finite q > 1 and any function b provided both b and b' are uniformly bounded. 
Now, as a consequence of the entropy balance (llObp , we have 


(152) 


(153) 


Divt 


g6s{gs,'ds) + Sgsiogi-ds), igsBigs^-ffs) + Sgslog{'ds))us +^Vooi^s 

Vs 

>0 in P'((0,r) X O), 


while (jl53p yields 


Curl; 


t^X 


b{'ds),0,0,0 bounded in L^((0,T) X O;R^^ ). 


Thus a direct application of the Div-Curl lemma 


gives the relation 


gs{g,'d)bi'd) = gsig,!)) 6(d), 


(154) 


for any 6 as in (11531) . 

Moreover, seeing that the sequence {e 5 s(p 5 , d 5 )d 5 }^^Q is bounded in Lp((0,T) x O) we deduce from 
(I155|) the formula 

gs{g,'&)'& = gs{g,'&)'&. (155) 

2. Parametrized measures and pointwise convergence of the temperature 

To show that p55l) implies strong convergence of the sequence we first remark that the 

approximate solutions solve the renormalized continuity equation 

9 t 6 (p 5 ) + div;j,( 6 (p 5 )u 5 ) + (^ 6 '(p 5 )p 5 - 6 (p 5 )^ div 3 ,(u 5 ) = 0 in T>'((0, T) x R^), (156) 

provided that gs and are extended by zero outside 17, and for any continuously differentiable function 
6 whose derivative vanishes for large arguments. 

Functions gs being square-integrable because of the artificial pressure in the energy equality, (11561) 
follows from p03l) via the technique of DiPerna and Lions [28] . 

Now it follows from (11561) that 

b{gs) —>■ b{g) in (^([O, T]; L(),g^j,(17)) for any finite p > 1 and bounded b. (157) 

Relations (I157|) and ()153l) yield __ 

g{g)h{d) = g{g) h{d), 

or, in terms of the corresponding parametrized measures (see Section |2?2T]) 

d^t,x{g,^) = M,x{g) ^ (158) 

Relation (I158|) means that oscillations in the sequences {£' 5 }, 5 >o {^< 5 }i 5 >o orthogonal, in the 

sense that the parametrized measure associated to can be written as a tensor product of the 

parametrized measures generated by {^> 5 } 5 >o- 

Let us consider a function H = H{t,x,r,z) defined for t S (0,T), a: G 17, and {r,z) G R^ through 
formula 

H{t, X, r, z) = r(^s{r, z) — s(r, 'd{t, x)j (^z — d^t, x)j . 
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Clearly, H is a. Caratheodory function: is continuous for a.e {t,x) € (0,T) x and 

H{-,-,r,z) is measurable for any (r, z) G 

Moreover, as both sr and sp are increasing functions of the absolute temperature, we have 


H{t, X, r, z) > x)J ^z — d(t, > 0. 


(159) 


f 

/O JQ. 


At this stage, we use Theorem 6.2 in [53], namely: weak limits of Caratheodory functions can be expressed 
in terms of the parametrized measure. Accordingly, we obtain 

lim / / (p{t,x)H{t,x, gs,'ds) dx dt = / / ip{t,x)[ / H{t,x, g,!}) dAt j;{Q,d)] dx dt 

Jq Jq Jq \J^2 ’ / 

ip(t,x)( [ gs{g,d){'d - d{t,x))dAt^x{Q,'d)) dx dt 

(p{t,x)( [ ps(p,d(t,x))(d - d(t,x))dAt,a;( 0 , d)) dx dt 
1 ^ ^ 

(p{t,x)( I gs(g,i9)(^ - i9(t,x))dAt,a:(g,^)) dx dt 
! ^ 7 r2 / 

= J j i^(^gs{g^ 'd)'d — gs{g, 'i9)i9'j dx dt = 0 for any (p G T’(( 0 , T) x fl, 

where we have used (I155I) to get the last equality together with (I158I) . observing that 

ps(p,'!?(t,x))(d - d(t,x)) dAt,x(g,^)= f gs{g,'d{t,x)) dAt^xig) f (d-■(?(<, x)) dAt,:j,(d) = 0 . 

Jr Jr 


f 

h JQ. 


f 

h Jn 


Jr^ Jr 

In particular, we deduce from (I159I) that -d^'d = which is equivalent to 

^ i? in L'‘(( 0 ,r) X n). 


(160) 


2.2.6 Pointwise convergence of densities 

1. The effective viscous pressure. Let us start with a result of P.-L. Lions m on the effective viscous 
pressure. 

In the present setting, it can be stated in terms of the parametrized measures as follows 

tpp{g, d)6(p) - ijjp{g, d) b{g) =TZ:[ip S]b{g) - TZ : [if S] b{g), (161) 

for any ip G 'DiTl), and any bounded continuous function 6 , where TZ = TZij is the Riesz operator defined 
by means of the Fourier transform: 


TZi 




— T~^ 


L|CP 






(162) 


Note that (11611) does not depend on the specific form of the constitutive relations for p, § and requires 
only validity of the momentum equation (11051) and the renormalized continuity equation (I156p for gg and 

U5. 


2. Commutator estimates. If viscosity coefficients v and g were constant, we would have 7^[§] = 
(|j 2 + 77 ) diva;(u), and (11611) would give immediately 



























where the quantity p— + ??)div 2 ;(u) is called the effective viscous pressure. 

In order to establish the same relation for variable viscosity coefficients, we write 

n : [V'S] = Ip diva;(u) + : [tpE>] - + 77 )diva;(u)|. 


Applying Lemma 4.2 in [?D], the quantity 


n : 


ip{v{'ds, Bi) (Va^u^) + r]{'ds, B5)diV:j,(u5)I^ 


- B^) + p{'ds, B5))diva;(u5), 


is bounded in the space L^( 0 , T ; for suitable 0 <w<l,p>lin terms of the bounds established 

in (11241) . (11251) and (11321) provided that the viscosity coefficients are globally Lipschitz functions of their 
arguments. 

If this is the case, one can use (11571) in order to deduce the desired relation 


: [ip S]b{g) -Ti:[ipS] b{g) = ip{-u + 77 )diva;(u) b{g) - ip{-v + ? 7 )diva;(u) b{g) = 


+ v)div^{u) b{g) - ip{-v + 77)div,j,(u) 6(p), 

where the last equality is a direct consequence of the pointwise convergence proved in (|146l) and (llBOp . 

If V and r] are only continuously differentiable as required by hypotheses of Theorem 12.11 one can 
write 

' B) = Y{d, B) + (1 - r(r?, B))i^(r?, B), 

(164) 

B) = Y{§, B) + (1 - B))r;(r5, B), 

where Y G T>(R^) is a suitable function. 

Now we have 

B5) (u^) + ry(d5, B5)diva;(u5)I = 

F(r?5,B5)z/(-d5,B5) (u^) + y(-d5, B5)7?(r?5, B5)diva;(u5)l+ 

|(1 - F(r?5,B5))z/('d5,B5) (u^) + (1 - ^(■da,B5))77(r?5,B5)div:j,(u5)l|, 

where the expression in the curl brackets can be made arbitrarily small in the norm of U‘{{Q,T) x 11), 
with a certain p > 1, by a suitable choice of Y (see estimates (11121) . (11191) . (11251) and formula (I133p . 
Thus relation (11631) holds for any v and 77 satisfying hypotheses of Theorem 12.11 


3. The oscillations defect measure. In order to describe oscillations in the sequence we use 

the renormalized continuity equation (1591) together with its counterpart resulting from letting <5 —^ 0 in 

(UMl). 

Although we have already shown that the limit quantities g and u satisfy the momentum equation 
(EB, the validity of its renormalization (1591) is not clear as the sequence {^?(5},5>o known to be 

uniformly square integrable and the Di Perna-Lions machinery |28j does not work a priori. 

In order to solve this problem, a concept of oscillations defect measure was introduced in [38]. 

To be more specific we set 


osc, 


([£<<5 e]{Q) = sup (limsup [ \Tk{gs) - Tk{g)\‘^ da; dt ), 

k>l ^ (5->-0 Jo ' 


(165) 


where Tk are the cut-off functions 

Tk{z) = sgn(z)min{|z|,/c}. 

As shown in [38], the limit functions g, u solve the renormalized equation (l59l) provided 
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• Qs, Us satisfy (11561) . 

• {'i<5}5>o bounded in L^(0, T; 




osCq[£i 5 —>• £i]((0,T) X fl) < oo for a certain q> 2. 


(166) 


Then, in order to establish (15^ it is clearly enough to show (11661) . 

4- Weighted estimates of the oscillations defect measure. In order to prove (11661) we use weighted 
estimates, where the corresponding weight function depends on the absolute temperature. 

Taking b — Tk in (116,11) we get __ 

p{g, 'd)Tk{g) - p{g, ■&) Tk{g) (167) 

= (^J^(i?,B) + 77(i?,B))div3;(u) Tk{g) - B) + 7?(il,B))diva,(u) Tk{g). 

As observed in (1451) . there is a positive constant Pc such that prig, "d) —Pcg^ is a non-decreasing function 
of g for any d. 

Accordingly we have 

p{g,d)Tk{g) - p{g,d) Tk{g) >Pc{g^Tk{g) - gi Tk{g)^. (168) 

Now let us choose a weight function w G C'^[0,oo) 

w(d) > 0 for -d > 0, w('d) = for -d > 1, (169) 

where a is the exponent appearing in hypothesis (1541) . 

Multiplying (jl67l) by w{d) and using (11681) we obtain 

w{f)[g^Tk{g) - g% (170) 

< + r]{d,B))diva:{u) Tk{g) - (^z/(-d, B) 77 ( 1 ?, B)]div 3 ,(u) Tk{g), 

Pc J Pc / 

where the right-hand side is a weak limit of 

Bs)+pids, B 5 ))div,(u 5 )) (Tkigs)-Tkig)) B)+ 7 ?(d, B))div,(u) (Tfc(p)-I^). 

Consequently, using the uniform weighted estimates (11191) together with hypothesis (IT^ and the growth 
restriction on w specified in (I169|) . we conclude that 


7(d) /4 


n L Pc 


(3'^(^>B)-H?7(-d,B))diva:(u) Tk{g) 


-^^^^(^i^{d,B)+p{d,B))divx{u)Tk{g) da: dt (171) 

Pc ^3 


< c li^inf \\Tk{gs) - 7fe(£')||L2((o_r)xn), with c independent of k. 

On the other hand, it was shown in Chapter 6 in [33] that the left-hand side of (I170p could be bounded 
below as 


[ [ w{d)(g3Tk{g) - gs Tk{g)) dx dt >limsnp [ [ w{d)\Tk{gs) - Tk{g)\‘^^ dx dt. (172) 

Jo Jn ^ ^ s^o Jo Jn 

Then (|170l) together with p71l) . p72L give the weighted estimate 

limsup / [ w{d)\Tk{gs)-Tk{g)\^ dx dt<c limmi\\Tk{gs)-Tk{g)\\L^((o,T)xn)- (173) 
<5->-o Jo Jn '*“*■0 
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Now taking q> 2 and we can use Holder’s inequality to obtain 



\Tk{Q 5 ) - Tk{Q)\‘‘ da; dt = 


\Tk{gs) - Tk{gW{l + dx dt 


(174) 


< c( [ [ iT’fcCe.s) - T’fc(£')|* (1 + i?) ^ dxdt+ [ [ (1 + i}) 2(8-39) da;dt'). 

^Jo Jn Jo Jn ' 

On the other hand, in accordance with (I112p . (11251) ■& S L°“(0, T; L'^(H)) O L^(0, T; L®(H)), and a simple 
interpolation argument yields ■& S L’'((0,T) x H), with r = ^. Consequently, one can find q> 2 such 
that 

39(1 -\- a) ^ 46 

2(8- 3g) 

provided a complies with hypothesis (l70)) . 

Then finally relations (11731) . (11741) give rise to the desired estimate 


limsup f j iu(i?)|rfc(£>5) - rfc(£<)|‘^ da; dt < c(^l + li^inf ||rfc(£i5) - Tfe(£i)||L 2 ((o,T)xn)), 
yielding (Ilfi6|) . 

5. Propagation of oscillations and strong convergence. According to (11661) we know that the renor¬ 
malized continuity equation (|59p is satisfied by the limit functions g and u. In particular 


dtLkig) + dYv,,{Lk{g)u) +Tk{g) diva;(u) = 0 in V'{{0,T) x R^), 

/* ^ T' (- 2 ^) 

provided g, u have been extended to be zero outside O, where Lk{g) = g — 5 — 

Ji z 

On the other hand, one can let d —0 in (11561) and one obtains 


dz . 


dtLkig) -I-divx(Tfe(e)u) + Tkig) diva,(u) = 0 in T>'((0,r) x 


(175) 


(176) 


Now, following step by step Chapter 3 in [42] we take the difference of (|175l) and (11761) and use (11671) and 
(I168P to deduce 


[ Lkig)iT) - Lkig)iT) dx < [ [ (diva;(u)) (Tfc(£>) - rfe(£i)') da; dt for any r G [0, T], 

Jn Jo Jn ^ J 

where the right-hand sides vanishes for fc —>■ 00 due to (I166p . 

Consequently 

6’log(£') = elog(e) in(0 ,r)xO, 

a relation equivalent to strong convergence of ’"^hich means in fact that 

in L^((0,T) X 0). 


(177) 


2.2.7 End of the proof 

Note that we have already shown that p, u satisfy the continuity equation (1571) as well as its renormalized 
version (l5^ . Moreover, it is easy to see that the “averaged” momentum equation (I147p coincides in fact 
with (l60l) in view of the strong convergence results established in (11461) . (11601) and (11771) . By the same 
token, one can pass to the limit in the energy equality (llOSp in order to obtain (|67)) . Note that the 
regularizing <5—dependent terms on the left-hand side disappear by virtue of the estimates (Blip . (|112p . 
and (11361) . 
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Similarly, one can handle Maxwell’s equation (11071) . Here, the only thing to observe is that the terms 


— --—curla;H 5 are bounded in L^((0,T) x fl), for a certain p > 1, 

(T{Qs, Vs, Ha) 

uniformly with respect to 6 (such a bound can be obtained exactly as in (11331) . 

To conclude, we have to deal with the entropy inequality (11061) . First of all, the extra terms on the 
left-hand side tend to zero because of (llll|) . (I112|) . (jl30|) and (jl51l) . Furthermore, it is standard to pass 
to the limit in the entropy production rate keeping the correct sense of the inequality as all terms are 
convex with respect to the spatial gradients of u, d and B. 

Finally, the limit in the logarithmic terms can be carried over thanks to estimate (11301) and the 
following result (see Lemma 5.4 in [2^ holds 


Lemma 2.2 Let H C be a bounded Lipschitz domain. Assume that 

ids ^ id in L^((0,T) x H) and log('d 5 ) —^ log{id) weakly in L^{{0,T) x H). 


Then id is positive a.e. on (0,T) x Ll and log(d) = log(d). 


3 Other models 

3.1 Polytropic models 

The magnetic Navier-Stokes-Poisson case where magnetic permeability /r is a pure positive constant was 
studied by Ducomet and Feireisl in |30) . 

In the absence of magnetic field a model of pressure law more adapted to cold plasmas of the type 

pig,id) = pe{g)+ idpi){g), (178) 

with associated internal energy e(p, d) such that 

ep = \{p-9pe), e^ = Cv{id), (179) 

P 

introduced by Feireisl [39] in the case of constant viscosities and generalized to temperature dependent vis¬ 
cosities and conductivity in [40] , has been extended to the MHD case (with Dirichlet boundary condition 
~ 0 constant resistivity) by Hu and Wang in [5T] : 

Theorem 3.1 Let H C be a bounded domain with boundary of class ^ > 0. 

Suppose that the pressure p = p{g, d) is given by (11781) . that Pe,P-d are functions of their arguments 
on [0, oo) such that 

Pe(0) = w(0) = 0, 

p'eip) > P^if) > 0 for all p> 0, 

Pe{p)<a 2 P^, p^{p) <a 3 {lp for all p>Q. 

Furthermore, we suppose that the magnetic permeability p and the electrical conductivity a are positive 
constants and that the transport coefficients 

K = k{i}), V = viP), — := A = A('d) := gidf) — cyiff), 

(J11 
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are continuously differentiable functions of their argument satisfying 



Ci(l+1?“) < K{d) < C2(l+1?“), 

(180) 


0 < X{d) < A, 

(181) 


0 <v< v(d) < V, 

(182) 

and 

0 < cv < cv{d) < cv. 

(183) 

for Cl, 

C 2 , A, ig, V, cv, cy > 0 and a > 2. 



Finally assume the homogeneous boundary conditions 


H|an = 0, u|an = 0, q- 11190 = 0 , 

and let initial data qq, (£iu)o, ’Oq, Bq be given so that 

go G {gu)o G do G L^(W^), Bq G (184) 

go >0, do > 0, (185) 

iQs)o = gosigo,do), — |(6iu)oP, {ge)o = goeigo.'&o) & (186) 

00 

and 

diva; Bo = 0 m X>'(11), Bq ■ nian = 0. (187) 

Then problem (I57)l - (IS51) possesses at least one variational solution g, u, d, B on an arbitrary time interval 
(0,T). 

This model, with boundary conditions 

B • n|ao = 0 , curlaB x = 0 , 
has been extended to the variable conductivity case 

0 < CT < cr(£i, d) < CT, 

for constants ct, ct > 0 and for a > 2, but under additional constraints on viscosities, by Fan and Yu [36]. 

3.2 Barotropic models 

When the temperature fluctuations can be neglected, one obtains the barotropic-MHD system for the 


density p, the velocity u and the magnetic field H in a bounded region ft 

dtg + diva,(pu) = 0, (188) 

dt{gu) + diva;(pu (g) u) = diva,T + j x B, (189) 

9tH + curia;(H X u) + curia,(Acurla;H) = 0, div^H = 0, (190) 

Ulan = 0, = 0. (191) 


where the pressure is p{g) = Ag'^ and the transport coefficients v, rj, X are positive constants. 
Hu and Wang proved the following existence result m 
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Theorem 3.2 Let il G be a bounded domain with boundary of class , <5 > 0 and 7 > 3/2. 
Assume that the initial data g^, (^?u)o and Bq be given so that 

go e L^{n), {gu)o G L\n-,R^), Hq G (192) 

go > 0, (193) 

and 

div.Ho = 0 m V\n), Holan = 0. (194) 

Assume finally that f,g € L°°(Q,T x LI). 

Then problem (11881) - (11911) possesses at least one variational solution g, u, H on an arbitrary time 
interval (0, T). 


Let us mention that the same result holds (R. Sart [86]) with different boundary conditions on H. 

Let us quote for completeness some other models recently introduced. 

Y. Amirat and K. Hamdache [1] studied a barotropic ferroffuid model with a magnetization dependent 
pressure p = p{g, M) including in the ffuid description an extra equation for the angular momentum. 

Finally X. Hu and D. Wang [53] considered density dependent viscosities, extending to MHD the 
result of D. Bresch and B. Desjardins m relying on analysis of the so-called BD-entropy. In this last 
case, presence of the radiative Stefan-Boltzmann contribution is no-more requested but additional 
constraints on the pressure and transport coefficients (including magnetic permeability) are needed. 


4 Singular limits 


As usual in hydrodynamical problems, when putting the system in non dimensional form, various non 
dimensional numbers [101] (Mach, Reynolds, Peclet etc...) appear, which may be small (or large) in 
certain regimes of motion. Studying the corresponding limits may simplify drastically the description. 
Typically low Mach number limits of compressible flows are expected to be close to their incompressible 
limit. In MHD flows, asymptotic regimes appear naturally when the so called Mach number, Alfven 
number, Froude number or Peclet number are small. 

The simplest case has been studied by P. Kukucka [65] , for which one expects convergence toward an 
incompressible limit. 

More precisely for a suitable scaling (see [65]) one considers the following non dimensional form of 
the MHD system (l9])- (fT4)l - ([T7|) - ([24|) where magnetic permeability p is supposed to be a positive constant 
(and we note A(p, ■!?,£!) := the magnetic diffusivity), to simplify the exposition. 


+ curl2,(B X u) -I- curia; (Acur^B) = 0, 


with 


5tp-|-diva;(pu) = 0, 

dt{gu) + diva;(pu 0 = diva;S -b x B -b 

d /Ma^ , Ma^ 1 , 

™ ' BP dx = 0. 


d r M a , 


dtigs) + diva;(£>su) 


Pe 


-div 


"'(I) 


= cr, 




d > ^ \ Ma § : Va;U + 
n 


Ma^ 

'AP 


IcurLHp - 


1 q-V.^ 

Pe d 


where <i> is an external potential (gravitation for example). 


(195) 

(196) 

(197) 

(198) 

(199) 

( 200 ) 
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In order to simplify the exposition, following [SS], we suppose that Mach number Ma and Alfven 
number Al are both small Ma = Al = e with 0 < £ <C 1 while Pe = 1. We also discard the external 


potential force ($ = 0 ). 

We assume boundary conditions 

B-n|aj^= 0 , Exn|as^ = 0, u- 11195 ^ = 0 , §n x = 0, = 0, ( 201 ) 

and initial conditions 

{B,g,u,'&){0,x) = {Bo,go,Uo,do){x) for any a: S (202) 

“ill-prepared” in the sense that 

Bo(a;) = B[|^](a;), goix) = g + e^'^lix), uo(a;) = u^^](x), -do(a;) = d + d[)^£(a;), (203) 

where 'g and p are two positive constants. 

Formally expanding (I195p - (jl99l) with respect to e we actually find at lowest order the incompressible 
MHD system with temperature 

dtB -f curl 2 ,(B X U) -f curia, (^( 0 , i?, B)curla:B) = 0, (204) 

div,(U) = 0, (205) 

g (^tU -f U • Va,U) -f Va,n = diva, (§(U, ^)) -b curla,B X B, (206) 

gcp{g,'d) ((9*0 -b U • Va;0) = diva, d,B) -b Va,0^ , (207) 

where cp{g^d) = e.^ + with boundary conditions 

® ~ 0 , U • njgQ = 0 , Va ;0 • = 0 , (208) 

and initial conditions 

B(0,x) = Bo(x), U(0, x) = Uo(x), 0(0, x) = i?o(a;) for any x e 11. (209) 

then one has the convergence result [65] 


Theorem 4.1 Let 11 C be a bounded domain with boundary of class , d > 0 and assume that 
properties of the state functions and transport coefficients given in Theorem \2.1\ are valid. 

Lei {pe, Ug, B£}e>o be a family of weak solutions to the scaled MHD system (11951) - (11991) ofi(0,T)xll 
supplemented with boundary conditions (12011) and initial conditions (12031) such that 

Qo,e weakly* in L°°{Ll), 

Uo,e —>■ weakly* in 

"do.e ^weakly* in L°°{n), 

Bo,e->B[|^^ weakly* in 

as £ —>■ 0. Then 

ess SMPt6(o,T)ll£>e(0 - < C'£, 

Ue^U weakly in 

B 

—-^B weakly in L'^(0,T;W^-'^(D;R^), 

£ 
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l9e-l9 


—>■ 0 weakly in i^(0, T ; R^), 


e 

where (U, 0,B) is a weak solution of the incompressible MHD system with temperature (12041) - (12071) with 
boundary conditions (12091) . initial temperature 


00 = 




cp{g,d) 






and initial velocity given by 


eU(0) = lim iJ[6ieUe](0), 


where H is Helmholtz’s projector on the space of solenoidal fields. 


The proof follows the stategy of [H] : 

1. One first proves uniform estimates for the sequence {pejU^, and the entropy rate 

2. One pass to the limit £ —?► 0 in the system, the most delicate term being the convective term 
QeUe (g) Ue, for which One uses the Helmholtz decomposition together with the associated acoustic 
equation. 

More complex situations have been investigated: A. Novotny, M. Rwzicka ang G. Thater m have 
considered the case Ma = e, Al = e, Fr = e, Pe = e^, Y.S. Kwon and K. Trivisa [55] have studied the 
case Ma = e, Al = y/e, Fr = ^/e, Pe = 1, E. Feireisl, A. Novotny and Y. Sun [43] have considered 
Ma = e, Al = Fr = Pe = 1 Re = l/e", for a n > 0 (large Reynolds number regime) and S. Jiang, Q. Ju, 
F. Li and Z. Xin [56] have studied more complex regimes. In all of these situations the target system 
is the incompressible MHD system. 

5 Strong solutions 

Previous Sections dealt with weak solutions, belonging to Sobolev spaces, for system ([3)-(0-([I7])-([2I]), 
together with ([7]) and ([ 8 ]) for the definition of the electric field and the current density. It is not known 
in general if these solutions are unique. 

Another approach is to look directly for strong solutions, which are assumed to be regular. In order 
to do so, the general method is to prove local existence of smooth solutions. This is done by applying 
general local existence results on hyperbolic and hyperbolic-parabolic systems (Subsection 15.11 below). 
Then, it is possible, assuming that the initial data are sufficiently close to an equilibrium state, to prove 
that the parabolic nature of the equation dominates. This implies a priori bounds which allow to prove 
that the solution remains in the neighbourhood of the equilibrium as long as it exists. This in turn 
allows to prove global existence. This method, due to Kawashima and his co-authors, is described in 
Subsection 15.21 Subsection 15.31 is devoted to existence and uniqueness results in critical (Besov) spaces. 

5.1 Local existence results 

In |95j , it is proved that a symmetric hyperbolic-parabolic system has a unique local solution. Therefore, 
its application to the present system only amounts to proving that the system is symmetrizable. The use 
of entropy variables is necessary to give a symmetric system for the fluid part (equations (im) - (IT71l - (I^ L 
See for instance m for the details. The part with the magnetic field is dealt with using the method of 
m, so the results of ^5] apply. The result may be summarized as follows: 
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Theorem 5.1 Let Q, = and assume that the thermodynamic functions 

p = p{g,'0), e = e{g,'d), s = sig,!}) 

are interrelated through (HSl), where p, e are continuously differentiable functions of positive arguments 
g, tt satisfying (l34)l - ((35ll - (IM)l - (l37)l . 

Furthermore, we suppose that each of the transport coefficients 

V = v(:d, |H|), p = p{d, |H|), kf = KF{Q,'d, |H|), a = a{g,d, |H|) and p = /i(|H|), 

are positive smooth functions of their arguments. Consider an equilibrium state u = 0, i?, B = 0 of the 
system (l^- (IT^ - (fT7ll - (l2T1l . and let the initial data data go, uq, ido, Bq be given so that 

{go-g,uo,do-^,Bo)GH^iR^). (210) 

Then there exists T > 0 such that the system (|n|) - (1141) - (1171) - (1211) . with the initial condition g^, )o, -do, 

Bq, has a unique strong solution on [0,T]. This solution satisfies g,T9 > 0 m x [0,T], and T depends 
only on (po “ Uq, i?o “ Bq) „3, inf o.nd inf po . 

" ^ ^ 11 -H R3 R3 

Some remarks are in order. 

First, the above result is stated in m, for electro-magneto-fluids (which, in addition to the above 
systems, includes equations for the charge Pe and the electric field E) with two-dimensional symmetries. 
However, the aim of m is to prove global existence, and this is why their study is restricted to two- 
dimensional flows. The local existence result is still valid for the 3D case. As a matter of fact, it is an 
application of the results of |95] . which are valid in any dimension. 

Second, the same result can be stated with instead of H^, for any s > 3 (see [23 Section 4]). 
Third, Theorem IS.ll is proved in (MJ Theorem 2.1] in the special case of constant coefficients p,kf,<J 
and p. The proof is not given, and is claimed to be an easy adaptation of that of [77]. Looking at the proof 
of [77], it is restricted to the case B = 0, that is, the compressible Navier-Stokes system with thermal 
conductivity. Their proof readily applies to the present case, and allows in fact for variable coefficients 
satisfying the above assumptions. A complete proof is also given in [37], in the case of a bounded domain 
(the method is easily adapted to D = K^), with possible occurrence of vacuum, that is, p > 0 instead of 

p > 0 . 

Note also that, in Theorem Ol we have assumed dMl), that is, the fluid is a perfect gas. Now, the 
result is still valid for a general equation of state satisfying (|35l)-(|36|)-((37|) only (see ED). 

Finally, the assumptions v,p, KF,cr, p > 0 may be weakened in the following way: it is in fact possible 
to treat cases in which some of these coefficients are identically 0. For instance, using the same method 
as in [3D, it is possible to deal with the case v, p,a, p > 0 , up = 0 , the case v, KF,cr, p > 0,p = 0, and the 
case v,a, p > 0 , Kp, p = 0. 


5.2 Global existence results for small data 

Theorem 15 .1 1 being proved, one can then turn to the proof of global existence of a smooth solution for the 
system under consideration. For this purpose, the key ingredient is an a priori estimate which we now 
detail. For this purpose, we define, for any 0 < < ^ 2 , the quantity 


Ns{ti,t2)= sup ||(p-p, U,1?- i?,B) (t)| 


i 1 ^t^t2 




r (iKVp, VB)f^._,(„3) (t) + ||(Vu, Vi?)||^.(R3) (t)) dt 

.J t\ 


1/2 


( 211 ) 


Then, differentiating the system if necessary, and multiplying it by some well chosen test functions, it is 
possible to prove the following kind of estimate: there exists 5 > 0 such that, if 


N,{Q,tf <C {N,{Q,Qf F N{Q,tf) , 
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for some constant C which depends on 6, but not on t > 0. This in turn aiiows to prove that there exists 
5 > 0 and C{S) > 0, both independent of t, such that 

^ 3 ( 0 , 0 ) <5 ^ N3{0,t) <C{S)N3{0,0). ( 212 ) 

Using (12121) . one proves a bound on 7 V 3 ( 0 ,t) for any time t of existence of the iocai in time solution. 
Applying the local existence result successively any each small time interval, one finally proves global 
existence: 

Theorem 5.2 Let U = and assume that the thermodynamic functions 

p = p{Q,-d), e = e{g,'0), s = s(p,d) 

are interrelated through (|25l). where p, e are continuously differentiable functions of positive arguments 
Q, § satisfying (|55|) - (I571) . 

Furthermore, we suppose that each of the transport coefficients u,r], KF,a, p are positive constants. 
Consider an equilibrium state p, u = 0 ,'i?, B = 0 of the system ([ ^ - (fT4l) - (fT7l) - (|2T|) . and let the initial 
data go, Ug, "do, Bq be given so that 

(po- 0 ,uo,do-d,Bo) (213) 

Then there exists 6 > 0 such that, */iV 3 ( 0 , 0) < 6 (N^ is defined by (12111) . then problem (l^- (IT4t - (|T7)) - (l2T1) . 
with the initial condition go, ug, dg, Bg, has a unique global strong solution (p, u,d,B) such that 


{g-g,B)inC° {R+,H^)nC^ , (u,d-d) G C° {R+,H^)r\C^ {R+,H^) , 

V (u,d-d) G . (214) 

This solution satisfies the estimate 

VtGK+, No(Q,t)<CNo{Q,Q), 


where C depends only on 6 > 0, and is defined by (12111) . 

This result is proved in [84l Theorem 2.5]. Considering the strategy of proof, it seems that the fact 
that the coefficients are constant is not important, although we did not check all the details. 

In [61) . the same result is proved for a more general system and with non-constant coefficients, but 
only with two-dimensional symmetries. Note that, due to the structure of their system, which includes 
an equation for the electric field and the density of charge, the equilibrium state does not necessarily 
satisfies B = 0. 

In essence, the above result is proved by studying the linearized equation near the equilibrium state 
(p, u = 0, d, B = O). This is a linear hyperbolic-parabolic system, for which energy-type estimates allow 
to prove the above result. Hence, a natural consequence is that the solution converges, as time tends to 
infinity, to the equilibrium state (p, u = 0, d, B = O). Such a convergence has been studied in [84], with 
the above assumptions. The convergence following convergence is proved in this paper: 


Vfc = 1,2,3, 


(e- 0 ,u,d-d,B)| 


L2(R3) < 


c 


(1 -(-1)2(9 2)+2 


(215) 


Vp > q, 


Here, q G [l, |) is such that (pg 


||(p-p,u,d-d,B)||^^ < 

^ ^ (l-bl)2U p) 

g, ug, 'do - d. Bo) G L‘>{R^). 


(216) 


As in the case of variational solutions, barotropic fluids have been extensively studied. In such a case, 
the system becomes (I188|) - (I189I) - (I190I) - (I191L with a pressure law p{g) = g'>', 7 > 1. All the coefficients 
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are supposed to be constant. In [71 [37], local existence is proved in the case A = 0, that is, when the 
parabolic term in (11901) vanishes. The main idea of the method is to fix the velocity field u, and then 
solve the problem with respect to all other variables. Then the Schauder fixed-point theorem is applied. 
Note that in [97], g !->■ p{q) is only assumed to be strictly convex and increasing. Global existence is 
proved in |14j under the assumption of small initial data, using a method similar to the one explained 
above. m proposes a generalization to the case of a system with Coulomb interaction, and [99] to the 
case of a bounded domain. The result in m is similar, but allows for vacuum in the system of equations. 
Finally, [H] gives a similar kind of result, but with an norm to estimate the initial data. Note finally 
that the asymptotic behavior as t —>■ -foo is also proved in [la 1111199] with the same kind of results as 
above. 

5.3 Critical regularity 

After the works of R. Danchin [22] [23] [24] on compressible Navier-Stokes system in a critical functional 
framework, some recent works have been devoted to the extension of his results to the compressible MHD 
system. 

For the reader’s convenience, let us first recall a few basic definitions concerning Besov spaces (a much 
more detailed presentation can be found in 1^). 

Let us first introduce the homogeneous Littlewood-Paley decomposition. Let y : K" —>■ [0,1] be a 
smooth nonincreasing radial function supported in i?(0,1) and such that y = 1 on B{0, 1/2), and let 

■= x(C/2) - y(0- 

Then we define the homogeneous Littlewood-Paley decomposition (Afe)fcgz over R" by 

AkU := (p{2-'^D)u = 


where J- stands for the Fourier transform on R". 

For any s G R and (p,r) G [l,-|-oo]^, the homogeneous Besov space i/p,.(R") stands for the set of 
tempered distributions / such that 




< oo. 


and we also note R® := i?| 2 - 

One has the formal decomposition 

u = ^ AkU. 
feez 

Owing to the fact that the solutions that we consider here do not have the same regularity in low and high 
frequencies, it is useful to use hybrid Besov spaces by splitting any distribution u into its low frequency 
part and its high frequency part [^ as follows: 

= and u^ = ^AfeU. (217) 

k<0 fc>0 


Then we define for s, t G R 


||/||b=. := ^ 2®'=|| A,/|U. + ^ 2‘'=|| A,/|U., 


and the hybrid Besov space B®’* 


k<0 k>0 

is defined for m = — [d/2 -|- 1 — s] by 
:= {/ G 5'(R®*) : ||/||b.,. < oo}. 
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if TO < 0 , or 


B*.* := {/ e S'(R‘^)/Vm : ||/||b^.‘ < oo}, 

if TO > 0, where Vm is the set of polynomials of degree less or equal to to (of course in any case := i?®’'*). 
Let us first rewrite the barotropic system 

dtg + div3;(pu) = 0, (218) 

dt{gu) + diva:(£<u 0 u) - VxV{g) = liAu + A)Va:div 3 ;(u) + H • Va;H - iv 2 ;(|Hp), (219) 

9tH + div^(u) H + u-V:,H) =^AH, div^^H = 0, (220) 

{g, u, H)t=o = [go, uo, Ho), diva,Ho = 0. (221) 

where the pressure is p[g) = Ag^ and the transport coefficients v,ri,X are positive constants. 

Hao proved the following existence result [15] 


Theorem 5.3 Let d > 3, p > 0, 2fj, + dX > 0, v > 0 and h an arbitrary constant vector. 

Assume that the initial data satisfy go —f G ^ ^ )o e g 

Then there exists a small number e > 0 and a constant M such that 

ll£’o ~ p||s<i/2-i.<i/2 + ||uo||5d/2-i + ||Ho||5d/2-i < e. 

then system I|218I1 - (I221|) has a unique global solution on an arbitrary time interval (0,T) such that [g— 
p, u, H - h) S E, with 

E := ^ (' 5 d/ 2 )d+d^ ^ ^d/ 2 +i.d/ 2+2 ^ ^^d/2+i^d+d^ ^ 

and satisfies 

||(p — P, U, H — h)||£; < M(||po — p||Bd/ 2 -i,d /2 + ||uo||Bd/ 2 -i + ||Ho||Bd/ 2 -i) . 

In the case of a polytropic fluid with constant viscosities the system rewrites 

dtg + div 3 ,(pu) = 0, (222) 

dt[gu) + div,„(pu 0 u) - VxP[g) = /^Au + (p + A)Vxdiva;(u) + H • Va,H - iVa;(|Hp) + pf, (223) 

i9t(pi?)+diva;(pui?)-KAr?+pdiva;(u) = ^|Va,u+V^u| = p(div,cu)^+cr(curla;H)^+H-Va;H-iVa;(|Hp)+pf, 

(224) 

atH + div^(u) H + u-V,,H) = i/AH, div^H = 0, (225) 

(p, u, H)t=o = [go, Uo, Ho), diva;Ho = 0. (226) 

Bian and Guo proved the following local existence result |7] (see also m 0) 


Theorem 5.4 Let d > 2, p > 0, 2p + dX > 0, v > 0. 

Assume that the initial data satisfy po — p S , Uo G [Bp[lY, i?o G B^^^ and Ho G [Bp^J’)'^, with 
go — p bounded away from 0 . 

Then there exists a positive number T such that the system (12221) - (12251) has a unique solution (p, u, rl, H 
on the time interval (0, T) such that 

p-pGC([ 0 ,r);<t), 
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u e (c{ [0, T); B^Jf- ') n ([0, T); +')) ^ , 

e C{[0,T);B:^^f-^)nL\[0,Ty,B^^/f), 

H G (^C{[0,T);B^Jf-^)nL\[0,T);B^Jf+^)y. 
Moreover if p G [2,d] the solution is unique. 


The first step of the proofs of both of these critical results is to use a smooth approximation of the 
data. For this regularized problem, previous results of R. Danchin [23] allow to prove local existence 
using a linearization argument. Then, uniform bounds are proved on the solution, allowing to pass to 
the limit as the regularization parameter vanishes. 


6 One-dimensional models 

In this section, we give an account of results in the one-dimensional case, which allows for drastic simpli¬ 
fication. In such a case, the use of Lagrangian variables gives a much simpler system which is amenable 
by standard techniques. 


In dimension one, it is commonly assumed that the unknowns depend only on the first coordinate 
of xi of X. The velocity of the fluid is still a three-dimensional vector, but its first component plays a 
special role. We denote it by u, while the transverse velocity is w: 



Hence, denoting B = ( 6 i, 62 , the constraint div(B) = 0 gives 9i6i = 1, so that hi is independent 
of xi- Further, if / is a vector field which is independent of cci, then the first component of curia,(/) 
vanishes. Hence, © implies that dtbi = 1. As a consequence, we may assume without loss of generality 
that hi = 1. Hence, denoting by b = (& 2 , ^ 3 )^ the transverse magnetic field, that is. 



we find the following system 

' dtQ + dxi (gu) = 0 , 
dt (gu) + dxi {gu^ +p) = 9a,i 


V + 'n] dxiU 


-h - dx 


\f |bP \ 1 


74 \ ' 


= dxt 

( 1 Mda,aU 


dt (pw) -I- dxi ( guw - - ] = dxi 


dt£ + dxi 


5tb -k dxi (uh - w) = dxi ^ , 


+ dxi {npdxt'd) 


(227) 
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where the total energy £ is given by 

^ = \s + |wp) + ge + M . 

Recall that is defined by m, and that, if /r is a constant, we have 


M 



2/x 


In what follows, we will always assume that /x is a constant. 

In the special case of a barotropic fluid, that is, without the energy equation (fourth equation of (12271) 
and with a pressure law given by p{g) = Ag^^ 7 > 1 , it is proved in [ 102 ] that the above system has 
a unique global smooth solution in the spatial domain ft = ( 0 , 1 ), for initial data (po, uo, wq, bg) such 
that go > 0, go G , (uq, wq) G Hq D bg G Hq. Further, a compatibility condition is needed, in 
the spirit of those used in the case of Navier-Stokes equations [T7J [THl EJ). All the coefficients v, g, a, fi 
are assumed to be constant. The method of proof is here again to prove local existence using Banach 
fixed-point theorem. Then a priori estimates are proved, which allow to apply local existence again and 
extend the time interval to infinity. These a priori estimates give, in a addition, a stability of the solution 
in Sobolev spaces. 

Next, in the paper [103] . Zhang and Xie study the full system (12271) . with additional gravitation terms: 


' dtg + dx^ (gu) = 0 , 
dt (gu) + [gv^ +p) = 


dt (gw) + ( guw - - ] = 


v + v] dx^u 


-h ■ dx^[^] + gdx^ip, 


|bP \ 1 


74 \ ' 


= dxi 

( + ?7 1 udx^u 


dt£ + dxi 


+ dxi • dxi + ^ 2:1 • dx^w) + gudxiiJ, 

dth + dxi (ub - w) = dxi ’ 


-b dxi iKFdxi'&) 


(228) 


In this system, ij) is the gravitation field, defined by Poisson’s equation: 


-9l,'tlJ = Gp inf]. 
Ip = 0 on d£l. 


(229) 


where G > 0 is the gravitational constant. The conductivity up = Kp{g,'d) is not assumed to be a 
constant, but satisfies the growth assumptions 


^(l + i?'^) < KF, 


dup 

dg 


< l + r?«. 



(230) 


for some constant G > 0. Further, the equation of state includes the radiative term -d'^, that is 


p = Rp'd+^'d^, 


e = Gvi?-b-z?^ 

g 


(231) 


The method is similar to [102] . This result has been improved by Qin and Hu [85] , in which q > 
0.3 is allowed. 
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Another approach in this one-dimensional setting is to use Lagrangian coordinates. Here, we assume 
that we use a bounded domain, which, without loss of generality, we may assume to be (0,1). We also 
assume that po(xi)dxi = 1. Setting 

y = y{t,xi)=[ (232) 

Jo 

we then have 0 < y < 1 , and, due to mass conservation 

/ g{t,xi)dxi = / go(t,xi)dxi = 1. 

Jo Jo 

Thus, defining v = Ijp the specific volume of the fluid, we immediately have, using Dirichlet Boundary 
conditions 

dty + gu = 0. 

Hence, system (12271) becomes 

" dtv — dyU = 0 , 


dtu + dyp = dy 


\v + P 


dyU 


-h- dy{ — 


5tW - I - I =d. 




dt {v£) + dy [pu+ - w • b ) = d. 


dt (r;b) - dyW = dy i —dy ( - 


+ rj] H- dy ( — 1 + -w • d, 

V V (TV \p, J V 




av \p. 


(233) 

In principle, this system is set on the interval H = (0,1) with homogeneous Dirichlet boundary conditions: 


(m, w, b, 9yi?) = 0 on (234) 

However, the corresponding Cauchy problem (H = R) is also studied by some authors, with the same 
kind of boundary conditions. System (I233I) - (I234I) is completed by initial conditions (wq, mq, wq, do, bo). 

The first existence result concerning the Cauchy problem (12331) is the global existence of Kawashima 
and Okada [62]. The equation of state is assumed to satisfy natural conditions similar to (l34l) . It is 
proved in [52] that, given a constant state (u, u = 0 , w = 0 ,d, b) , (with u > 0 and d > 0 ) if the initial 
condition is close to this state in 7J^(R), then system (I233|) - (I234I) has a unique global smooth solution 
{v, u, w, d, b) such that, for any t > 0, {v{t) — v, u{t), w(t) — W, d(<) — d, h(t) — b) S i7^(R). Here again, 
local existence is proved using classical results, and global energy estimates are proved for these smooth 
solutions. This allows to extend the existence result to any positive time. The result may be summarized 
as follows: 


Theorem 6.1 (Kawashima-Okada, 1982 [62]) Consider system (12331) set m H = R. Assume that 
the equation of state satisfies 


dp 

dg 


> 0 , 


de 


> 0 , 


and that de = ddS — pd j , where S is the entropy of the fluid. Assume that p is a positive constant, 
and that p, v, up are positive smooth functions of g and d. Let 


(u, u, w, d, b) = (n, 0,0, d, O) 
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be a constant equilibrium state of (|233D . Let (tip,M q, wq, bg) be an initial data such that 

{vo - V, uo, Wo, -do - -d, bo) G 


If ll'i'o “ v,uo,wo,'do — i?,bo||^ 2 (R) *5 sufficiently small, then (12331) has a unique global smooth solution. 


Actually, coefficients are allowed to vanish, provided they are constant: the situation allowed are (a) 
ly = rj = 0, yL, kf > 0, (h) iy,r], yL > 0, Kp = 0, (c) rj = v = Kp = 0, fj, > 0. 

Another result can be found in m- In this contribution, system (12271) is considered on the domain 
Ll — Ll{t) = {0,x{t)), where the free boundary x{t) is defined by x'{t) = u{t,x{t)). The boundary 
conditions are given by 


P - [ + da:^U 


= 1 . 


(235) 


(w,b,a2,id)afj(t) = 0, M(t,0)=0, 

Xi—x{t) 

Then, the Lagrangian coordinates are used to transform the system into (12331) . in the domain LI = (0,1), 
with the boundary conditions 


(w, b, = 0 , ■u(t, 0 ) = 0 . 


+ dyu 


= 1 . 


J y=l 


(236) 


We are back to a system which is similar to (j233l) - (l234l) . to which the strategy of Kawashima-Okada 
[52] is applied, thereby proving global existence if the initial data is close to a given constant state, in 
Sobolev spaces. Moreover, in the special case of a perfect gas equation of state, a global existence of weak 
solutions is proved. For this purpose, the initial data, which may be discontinuous, is approximated by 
smooth functions. This allows to apply the global smooth existence result. The estimates established 
in the smooth case are then used again to prove that one can pass to the limit as the regularization 
parameter tends to 0. This is based on Aubin-Lions Lemma (see min]). All these results are proved 
using the assumption that a = 1 , 77 , zz are bounded smooth functions of v, that the energy and pressure 
of the fluid satisfy some polynomial growth property as functions of d, and the fact that 


KF{v,d) + |9„kf(u,7?)| + \dlKF{v,'d) \ < C (1 + i?'^), 

for some constant C > 0 and some power q > 2. This is in sharp contrast with [62] . in which up is 
assumed to be constant. 

Using a similar method, Wang [96] proved that system (I233I) - (I234I) has a global weak solution in 
H^(L1). The same kind of assumptions are used on the equation of state and the diffusion coefficient Kp. 
In a similar setting, it is proved in [16] that the solutions depend continuously on the initial data. 

In the case of constant coefficients, with an equation of state satisfying 


dp, . 


3 ( 7 ?) ■ 




< c, 


the Cauchy problem (]233D set on fl = M is studied by Hoff and Tsyganov in [351 EO]- It is proved that 
weak solutions exist and are unique, provided the total energy is small. Under an additional integrability 
assumption on dth, continuity with respect to the initial data is proved (in Lf{[Q,T] x R), for any T > 0). 

In [35], for system (I233I) - (I234I) . the limit of vanishing shear viscosity —>■ 0 is studied, under the same 
assumptions as in m- It is proved that the solution converges to that of the limit problem strongly in 
L^([0,T] X LI) for p, w, 7 ?, and strongly in L'^{[0,T], H^{Ll)) forw,b. 

Finally, still in Lagrangian coordinates, but in the special case w = 0, system (I233I) - (I234I) has been 
studied by Iskenderova [5S]. Assuming constant coefficients and a perfect gas equation of state, Isk- 
enderova proves that the system has a unique weak solution. Note that the initial density po is allowed 
to vanish, provided that dy {po^'^ £ and that 

Po{,y)<C jni pq{z), 

0<z<y 
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for some constant C > 0. Here again, the method of proof is a regularization of the system: existence 
and uniqueness is easily proved for a system in which po ^ > 0- Then, the limit e —>■ 0 is studied using 

a priori estimates. 
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